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ABSTRACT 

Using three-dimensional, moving-mesh simulations, we investigate the future evolution of the re- 
cently discovered gas cloud G2 traveling through the galactic center. We consider the case of a spherical 
cloud initially in pressure equilibrium with the background. Our suite of simulations explores the fol- 
lowing parameters: the equation of state, radial profiles of the background gas, and start times for 
the evolution. Our primary focus is on how the fate of this cloud will affect the future activity of Sgr 
A*. From our simulations we expect an average feeding rate in the range of 5 — 19 x lO^^M© yr~^ 
beginning in 2013 and lasting for at least 7 years (our simulations stop in year 2020). The accretion 
varies by less than a factor of three on timescales < 1 month, and shows no more than a factor of 10 
difference between the maximum and minimum observed rates within any given model. These rates 
are comparable to the current estimated accretion rate in the immediate vicinity of Sgr A*, although 
they represent only a small (< 5%) increase over the current expected feeding rate at the effective 
inner boundary of our simulations (r = Ib^Rs ~ 10^^ cm), where Rs is the Schwarzschild radius of 
the black hole. Therefore, the break up of cloud G2 may have only a minimal effect on the brightness 
and variability of Sgr A* over the next decade. This is because current models of the galactic center 
predict that most of the gas will be caught up in outflows. However, if the accreted G2 material can 
remain cold, it may not mix well with the hot, diffuse background gas, and instead accrete efficiently 
onto Sgr A*. Further observations of G2 will give us an unprecedented opportunity to test this idea. 
The break up of the cloud itself may also be observable. By tracking the amount of cloud energy that 
is dissipated during our simulations, we are able to get a rough estimate of the luminosity associated 
with its tidal disruption; we find values of a few 10^^ erg s~^. 

Subject headings: galaxies: active — galaxies: ISM — Galaxy: center — Galaxy: nucleus 



1. INTRODUCTION 



Sgr A*, the supermassive black hole at the center of the Milky Way galaxy, is the most underluminous black hole 
observed. It has a bolometric luminosity of < 10^^ erg s~^ ( Narayan et al.|[l998 ), whi ch is less than 2 x 10 ~^ of the 
Eddington luminosity, LEdd = 5.2 x 10"^"^ erg s-\ for a Mbh = 4.3 x 10" black hole ( [Gillessen et"aLl|2009[ ). Part of 
the reason it is so underluminous is that it is starved for material. At present, its accretion is led by h ot [ksT ^ 1.3 
keV) gas, p ossibly expelled in the winds of luminous young stars that are orbiting within its vicinity ( Krabbe et al. 
[l991; Melia| |l992| j Baganoff et al."200 3j. The elec tron number density of this gas. He ~ 26 cm ^, suggests a Bondi 



deeding rate oT^nr^M© yr"^ ( ^Yuan^et al.||2003 
of polarization, along with F^r^d^y rotat ion ^rgi 
range 2 x 10"^ < M < 2 x IO'^Mq yr"^ (lAitken et 



at TBondi = GMbh/c^,6 



3.4 X 10^^ cm. However, measurements 



ments^ constrain the accretiori rate at the black 



al.||20QQ 



Bower et al. 



2003 



Mar rone et al. 



ole to be in the 

2007[ ) , suggesting it is 

difficult to actually get the gas down to the black hole. 'I'his is consistent with a number ol recent studies following 
the gas from these stellar winds, which have identified a number of physical processes that can significantly reduce the 
accretion effi ciency, including expulsion due to hydrodynamic winds ( Quataert 2004 ) , hang up due to residual an^lar 
momentum ( Moscibrodzka et al. 2006), and heating due to electron thermal conduction (Shcherbakov fc Baganoff 



20101). Furthermore, some ol the stars may produce winds that are too fast to be captured efficiently (Cuadra et al. 
?008). 

Within this environment, Sgr A* may occasionally enjoy a relative feeding frenzy. Tidal disruption of stars passing 
close to the black hole is one way the accretion rate could be temporarily boosted. Such an event may have been 
recently observed in a z ^ 0.35 g^alaxy as the e xtremely luminous Swift transient event Sw 1644+57 (Bloom et al.| 



2011 Burrows et al.|[2QTT Zauderer et al.| 2011). Within the Milky Way, the star S2 has made the closest observed 
approach to Sgr A"^, passmg within 120 AU at its orbital pericenter ( Gillessen et al.|[2QQ9 ). Even so, to be disrupted. 
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a star would nee d to pass within ^ 1 AU of Sgr A*. Such close passes are only expected to occur approximately once 
every 10^ yr (see Alexander||2QQ5 , for a review). 

A second opt ion, and the on e that we explore in this paper, is tidal stripping of gas clouds passing through the 
galactic center. |Gillessen et al. (2012) recently discovered a^as cloud (d ubbed G2) of ^ 3 Earth masses on a highly 
eccentric orbit (e = 0.9384) approaching Sgr A* ( Schartmann et al.|2012 ) . The cloud is expected to reach pericenter in 
2013.5, at a distance of approximately 270 AU (c^peri = 4.0 x 10^^ era). Because the self-gravity of the cloud is negligible 



(Burkert et al. 



2012), it is certain to be tidally stretched. The main questions are: 1) How much will hydrodynamic 



interactions with the background gas disrupt the cloud; 2) How much of the cloud material will be captured by Sgr 
A*; and 3) On what timescale will that captured material accrete? 
To answer these questions we have performed three-dimensional hydro dynamic simulatio ns of the cloud's infall. 



Simulations of this cloud have been done previously (Burkert et al.] |2012 Schartmann et al. 2012]), but only in two 



dimensions. Two-dimensional simulations offer the advantage that they allow one to cover the enormous dynamical 
range between the size of the cloud and the size of its orbit at reasonable computational expense. However, such two- 
dimensional simulations are highly artificial. One reason is that, in a Cartesian grid, a two-dimensional "cloud" has 
an infinite extent in its third dimension. Furthermore, the enforced symmetry ignores the effects of tidal compression 
along this third dimension. Finally, hydrodynamic instabilities, su ch as Kelvin- Helmho It z and Ray leig!:h- Taylo r, are 



Young et al. 12001 ) 



known to behave differently in two dimensions than in three (e.g. 'Sharp 1984 



Fritts et al. 



1996 



Kane et al. 



2000 



We were able to get around the problem of the e xtreme dynamical range by using the moving mesh capability of 
Cosmos ++ ( Anninos et al.||2005 Fragile et al.||2012 ). This allows our computational grid to follow the cloud along its 
orbit. Because we only need to resolve the cloud and not the volume of its entire orbit, we are able to proceed in three 
dimensions at reasonable computational expense. 

The paper is organized as follows: In Section [2j we review the most pertinent details about the cloud's structure, 
internal properties, orbital parameters, and environment. These are necessary for setting up the initial numerical 
problem, and also for estimating the relevant timescales for this problem, which we do in Section [3j This provides an 
opportunity to discuss the physical processes that will play a role in the disruption of the cloud as it passes pericenter. 
As it turns out, many of the timescales for these processes are shorter than or comparable to the orbital time. The 
competition among physical processes illustrates why a numerical approach is essential. In Section |4j we present the 
results of our three-dimensional, moving-mesh simulations following the cloud G2 through pericenter passage. We 
summarize our findings in Section 5] Additionally, since this is our first introduction of moving meshes in Cosmos++, 
we provide a brief description of the equations in Appendix [A] 

2. CLOUD PROPERTIES 

[Gil lessen e t al.| ( 1201 2*) detected the dense, dusty, and ionized gas clump, G2, with a dust temperature of 550 K and 
a gas temperature of order 10^ K, embedded in the diffuse 10^ K gas of the galactic center. They adopt an effective, 
spherical cloud radius in 2011.3 of 15 mas, which corresponds to Rc = 1.875 x 10^^ cm. The observed Br7 luminosity 
implies a cloud density of pc ^ 6.1 x 10~^^/y g cm~^, with a corresponding cloud mass of Mc « 1.7 x 10^^/y 
g, or approximately 3 Earth masses, where /y < 1 is the volume filling factor. We will use this simple spherical 
approximation along with a filling factor /y = 1 and a uniform cloud density to initialize our simulations, even though 
observations su ggest that G2 might, in fact, be the compact head of a larger, more diffuse distribution of cold gas 



( Burkert et al.|[2012[ ). Since the cooling timescale of the cloud is much shorter than its orbital period (see Section 3|, 
cooling plus ionization by the strong UV field of the central cluster of massive stars should combine to maintain the 
cloud temperature at Tc = 10^ K, which we take as its initial value, although we also explore two non-isothermal models 
that are allowed to evolve away from this temperature. Along with the ideal gas law, these parameters completely 
specify the cloud properties. 

We consider two m odels for the backgrou nd me dium. The first is based on the background parameters used in 
Burkert et aT] ( |2012[ ) and , Schartmann et aL] ( |2012| ), i.e. 



Pb = 1.3 X IQ-'^^T] 



Tb = 1.2 X 10^ 



<^1995 



1995 



g cm 



K. 



(1) 



(2) 



where 6^1995 = 10^ i?^ was the approximate distance of the cloud from Sgr A* in 1995.5 and Rs = 2GMbh/c^ = 
1.3 X 10^^ cm is the Schwarzschild radius of the black hole. We have assumed a mean molecular weight of ja = 0.614. 
The parameter r] < 1 is meant to acco unt for any unresolved X-ray sources contributing to the observed background 
X-ray luminosity (Sazonov et al.||2012 ). In practice, we use this coefficient to adjust the background density to ensure 
that the background is initially m pressure equilibrium with the cloud. 

One problem with this background is that it is convectively unstable, as the entropy, S = T j p^l'^ ^ decreases radially 
outward. In order to prevent the cloud fr om b eing disrupted by this co nvectively unstable atmosphere, we follow the 



same procedure that Burkert et al. ( 2012| and Schartmann et al. (2012) used to artificially stabilize it. This starts by 
additionally evolving a passive tracer field T. Fhe tracer is initialized to the same value as the cloud density inside 
the cloud and zero in the background. This allows us to later distinguish between zones with pure atmospheric gas 
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{T / p < 10~^) from those with significant cloud material {T/p > 10~^). Also, because the units we adopt in our 
calculations normalize the total cloud mass to unity, the value of the tracer field is a measure of the zonal mass relative 
to the total initial cloud mass. At the end of each computational cycle, zones identified as pure atmosphere have their 
densities and temperatures reset to the background values given in equations ([T]) and ^ and their velocities reset to 
zero. We tested other values (10~^ and 10~^) for the tracer floor and found that the mass accreted onto the black 
hole increases only by about 3% with each factor of 10 that the tracer floor is lowered. 
The second bac kground we consider is based on a more careful fit to the galactic center environment profile from 



Yuan et al.| (|2003|), designed to match Chandra X-ray observations. This gives a background of the form 

/ ^ X -1.125 

Pb = 1.3 X 10-^^7/ I ) g cm-^ (3) 



n995 , 
/ J \ -0.75 

T, = 10M-^) K. (4) 

\«1995 / 

Unlike the atmosphere in equations ([ij and this atmosphere is convectively stable. Otherwise this background is 
implemented in the same way as the nrst one. 

In both backgrounds we assume hydrostatic equilibrium within the gravitational potential of Sgr A*. This is certainly 
an oversimplification. Many physical processes are likely to play a role in determining the detailed structure of the 
galactic center environment. However, as there is no detailed data available at this time to more faithfully model this, 
we focus only on the dynamical interaction between the cloud and these simple atmospheres. 

We assume that the cloud is following a Keplerian orbit. We start the cloud either fr om 1995.5, shortly before its 
first detection, or 1944.6, the date corresponding to the apocenter of its current orbit (Burkert et al.||2Q12 ). In both 
cases the cloud is assumed to initially be in pressure equilibrium with the background. 

3. TIMESCALES 

In this section we consider the timescales associated with various physical processes that may act within the cloud 
or between the cloud and the background. We do this to estimate their relative importance to the evolution of G2. 
All estimates are scaled to the parameter values we use at the start of most of our simulations. This corresponds to 
the cloud position in 1995.5, although we use the size and mass of G2 as measured in 2011.3. This is appropriate since 
we initialize all of our simulations using the 2011.3 cloud properties. Thus, the timescales we calculate give us a good 
estimate of the relative importance of each of these processes near the start of our typical simulations. 

First we consider radiative cooling. The associated timescale is 

An2 

— »'Ay'(oi3)(i^)(^)"'(|)"'- (^) 

where the numerator is the energy density of the cloud (assuming a gamma-law equation of state with adiabatic index 
r = 5/3), the denominator is the volume emissivity, ric = PcKp'^h) is the number density of the cloud, ttih is the 
mass of hydrogen, /c^ is Boltzmann's constant, and p* = 6.1 x 10~^ g cm~^ is the initial cloud density used in our 
simu lations. To get a gene ral estimate of the cooling timescale, we use the same cooling rate A* = 3 x 10~^^ erg cm^ s~^ 



as in jBurkert et al.| ( [2012D . Since our estimated cooling timescale is much shorter than the orbital period of Torb 

= 138 

yr, we are justified in assuming that the cloud will maintain its initial temperature of Tc = 10^ K. 

Small gas clouds like G2, embedded in a hot environment like the galactic center, will lose mass due to evaporation 
as a result of thermal conduction. Whenever radiation and magnetic field effects are unimportant, the ev aporation 



timescale, in the case when the heat flux reaches its limiting value (saturated state), can be estimated from ( |Cowie fc 
McKee|[T977| ) 



-^^evap 



Pc\ f Pb\ ^ f Rc\ f \ 



where 

*4>,.1.7xlO.»(|)(|)'(j^^)"%s-. (7) 

is the saturated mass loss rate due to evaporation, pi = 1.3 x 10~^^ g cm~^ is the background density at digg^ = lO^i^^ 
for our default background, and = 1.875 x 10^^ cm is the initial radius of our simulated cloud. This evaporation 
timescale is somewhat shorter than the orbital period. It is, therefore, unclear if 02 could traverse the full distance 
from its apocenter. Further, it seems clear that once 02 breaks up as it approaches pericenter, the resulting fragments 
will likely dissolve quickly. 



4 



G2 will also experience ram pressure as it traverses the galactic ce nter. Ram pressure w ill compress the cloud on 
its leading edge and slow its motion relative to the background gas (Murray & Lin 2004 ). The importance of ram 
pressure, compared to the thermal pressure of the surrounding gas, is given by the ratio '^^ei/^s 6' where 



Cs,b = 



1/2 



=i.3xio^f^^y 

V0.614y 

is the sound speed of the background gas and 



■1/2 



1.2 X 10^ K 



1/2 



cm s 



(8) 



2GMbu 



:2.8 X 10^ 



df 



1/2 



Ml 



BH 



4.3 X lO^M,;, 



1/2 



df 



1995 



^apo 



1995 



1/2 



cm s 



(9) 



is the relative velocity between the cloud and the background. For simplicity, we utilize the free-fall formula. Noting 
that our background is stationary, we have 'Urei = '^c- Close to apocenter ((iapo = 1.3 x 10^^ cm), where 'Urei = '^0^0^ 
we see that gas pressure dominates and ram pressure is unimportant, while at the cloud's 1995.5 location, ram pressure 
begins to dominate. 

Since the self- gravity of G2 is not important, the cloud will also be subject to tidal forces from the background 
potential. Tidal forces will stretch the cloud in the direction of the black hole and squeeze it perpendicular to this 
direction. The net result is to redistribute the cloud along its orbital path. The acceleration of this stretching is given 

^RcGMbu 

ats . (10) 



This gives a timescale for tidal stretching of 



f2Rc 



:2.0 



1/2 



Ml 



BH 



4.3 X lO^Mc;^ 



-1/2 



1995 



3/2 



yr. 



(11) 



From this we expect the cloud to quickly be significantly tidally distorted. 

At the surface of the cloud, we have two different density fiuids moving relative to one another. This makes the 
surface susceptible to the Kelvin-Helmholtz instability. This instability will disrupt the surface of the cloud on a 
timescale 



rKU = {kVre\) 



1/2 



^9.0 



^rel 



2.8 X 10^ cm s-i 



1/2 



-1/2 



(12) 



where k is the wavenumber of the perturbation. We assume k ^ 1/ Rc to he the most disruptive wavenumber. This 
estimate suggests Kelvin-Helmholtz disruption may become important, especially as the cloud approaches pericenter. 

As the high-density cloud is ram-pressure decelerated by the low-density background environment, its leading edge 
will be subj ect to disruption by the Ray leigh- Taylor instability. In the linear regime, this instability grows on a 
timescale of ( Murray et al.|p!993 ) 



trt ' 



Mr 



irR'^kpijv'^ 



'11 



rel 
Mr. 



Pb^ Pc 

Pc - Pb 
1/2 



1/2 



1.7 X 1028 g 



-1/2 



-1/2 



'^rel 



2.8 X 108 cm s-i 



yr. 



(13) 



assuming the limit Pc^ Pb- Again the most destructive wavelength is k ^ 1/Rc. This timescale is comparable to the 
Kelvin-Helmholtz timescale indicating both processes may act on the cloud simultaneously. 
In summary, we find 

Tcool < rts < TKH < Trt < Tgvap ^ ^orb • (14) 
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The wide array of dynamical processes and timescales illustrate why realistic numerical simulations are so crucial. The 
only process that may be important that we do not capture in this work is thermal conduction, which in any case is 
not expected to be a factor until late times. 

4. SIMULATIONS 

The simulations use a Cartesian grid with a starting size of 14i?* on each side, resolved with 256 zones in each 
dimension, for a total of 256^ zones. Doing so, we achieve an initial linear resolution of Ax^y^z = 1.0 x 10^^ cm, 
comparable to that used in the 2D simulations of Schartmann et al. (2012). The fact that we can achieve a resolution 
in 3D comparable to that of a 2D simulation while covering the same effective volume illustrates the tremendous 
advantage of using a moving mesh. 

For a relatively simple case like the orbital motion of a gas cloud, it is most convenient to fix the motion of the mesh 
using the Keplerian velocity, rather than try to have the mesh dynamically adjust to the cloud's motion. Thus, we 
set the grid velocity (as defined in Appendix [A]) to Vg'y{t) = C'^c Kep(^)' where Vc^Kep is the Keplerian velocity vector 
of the cloud (assumed to lie in the x-y planej and ^ 1 is a parameter that allows us to account for the moderate 
slowing of the cloud due to ram-pressure deceleration. We find ( = 0.98 works well here. In the z direction, the mesh 
is given a small velocity toward the z = plane with a magnitude \Vg\ = 5i?*/(tperi — ^start), where tperi = 2013.5 and 
^start = 1995.5 or 1944.6. This compensates for the vertical slimming of the cloud due to tidal stretching and allows us 
to maintain an approximately constant resolution across the cloud in the z direction. This vertical grid compression 
stops at tperi- As the cloud passes pericenter, the grid motion in the x-y plane also changes, as we now wish to keep 
both the cloud and the black hole on the grid for the remaining evolution. To accomplish this, the x and y components 
of the grid velocity change to the following forms at t = 2014.7 and t = 2014, respectively: 

V,''=<K.,{t) /—_-f (15) 

i/max i/min 

where x^ax and 7/max are the maximum x and y grid positions during the time this grid velocity is implemented. This 
velocity field holds the "upper right" edge of the grid fixed, retaining the black hole on the grid, while allowing the 
"lower left" edge to stretch to follow the cloud. During this phase, the linear resolution of the grid decreases. We 
terminate the simulations at tgtop = 2020, with a linear resolution of Ax = 2.2 x 10^^ cm and Ay = 1.5 x 10^^ cm. 
The final z resolution is Az = 6.3 x 10^^ cm. 

The boundary conditions on the edges of the grid are set to outfiow conditions (all ghost zone quantities are set equal 
to the values of their nearest internal-zone neighbor, except that the velocity c omponent normal to the boundary is 
set to zero if it points onto the grid). Similar to |Burkert et al.[ (|2012) and Schar tmann et al.| (|2012| ), we also introduce 
a spherical "accretion radius" centered on the black hole, wnich is located at the origin of our coordinate system. 
All gas fiowing into this region is treated as being accreted. Each cycle, the amount of tracer material entering this 
region is tallied, and then all hydrodynamic variables inside it are reset to the static atmo sphere solution. For all o f 



our reported simulations, the accretion radius is set to race = 7bORsj half the value used in Schartmann et al. (2012), 
although we saw only modest differences (19%) in our measured mass accretion rate when we used the larger value. 

For this problem, we ignore the effects of the cloud and the background gas on the gravitational pot ential, and 
instead s i mply consider the (fixed) potential of the black hole. For this work we use the pseudo-Newtonian, |Paczyhsky| 



fc Wiita| ( |1980D potential ^ . . 

<'■' 

although at the distances we are considering, there is no practical difference between this and a standard Newtonian 
potential for a point mass. 

We experiment with three different equations of state: isothermal, isentropic, and polytropic. The first two are 
special cases of the polytropic equation of state, which specifies the pressure in the form P = tzp^ ^ where k is the 
polytropic constant, p is the gas density, and F is the adiabatic index. Our default EOS is the isothermal option with 
F = 1. This EOS, together with the perfect gas law (P = kspT/ prriH)^ ensures the cloud remains at a constant 
temperature. This can be thought of as representing a perfect balance between the heating of the cloud by the UV 
background and the radiative cooling of the cloud itself. For the isentropic case, we set F = 5/3, which allows the 
temperature to vary, but is held fixed. 

We also consider a third, more generic, polytropic gas with F = 5/3, but with allowed to v ary. This introduces an 



additional degree of freedom to the gas equations in the form of internal energy [equation (AS)]. The manner in which 
the pressure is recovered also differs in this case from the previous two options. For an isentropic (and isothermal) 
gas, the pressure is recovered directly from P = np^ ^ and the internal energy of the gas is not needed. By contrast, 
the polytropic option derives its pressure from P = (F — l)e, where e is the independently-evolved internal energy. 
For purely adiabatic fiows, polytropic and isentropic equations of state with the same n and F should give nearly 
identical results, as the pressure will change in time according to the adiabatic scaling. If there are shocks in the 
system, however, entropy is not conserved and the two equations of state may produce significantly different solutions. 

The equations of state are applied only to the cloud. The thermodynamic state of the background gas is specified 
by one of the two hydrostatic solutions in Section |2j together with the ideal gas law for the pressure. So long as 
cloud and background gases do not mix, one can apply the two procedures unambiguously. However, in reality cloud 
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Table 1 

Model Parameters 





Model 


EOS 


r 


Bkgd 


^start 


(10-8 Mo yr-l) 


cc 


_il_bl_95 


isothermal 


1 


1 


1995.5 


5.18 ±2.25 


cc 


_i53_bl_95 


isentropic 


5/3 


1 


1995.5 


7.75 ±6.15 


cc 


.p53_bl_95 


polytropic 


5/3 


1 


1995.5 


5.37 ±4.08 


cc 


_il_b2_95 


isothermal 


1 


2 


1995.5 


4.64 ±2.19 


cc 


.il_bl_44 


isothermal 


1 


1 


1944.6 


19.0 ± 12.1 



^ Mass accretion rate, time averaged from when accretion started until 
tstop = 2020. 



and background materials do mix as the cloud moves through the mesh, so a prescription for calculating a reasonable 
thermodynamic state for regions containing both gas and cloud materials must be provided. For our calculations we 
have chosen to impose a relaxation scheme that sets the temperature of mixed zones to a mass weighted combination 
of background and cloud temperatures. In practice this comes to: T^ix = {T / p)Tc + (1 — T/ p)Tg^ where Tc is the 
cloud temperature calculated from the appropriate equation of state, Tg is the background temperature calculated 
based on the zone position relative to galactic center, and T / p is the mass fraction of cloud material in the zone. We 
define mixed zones as those regions where 10~^ < T/ p < 0.8. Once the temperature is calculated in this fashion, the 
pressure follows from the ideal gas law. 

Table [1] summarizes the five models we co nside r in this paper, exploring three main parameters: equation of state, 
b ackground model, and star t date. In Section [TT] we describe our generic results and compare with the 2D simulat ions 



of Schartmann et al. 



2012). In Section 4.2 we compare results from the different equations of state. In Section 



4.3 



we explore our two different backgrounds. In Section [441 we consider the two different start dates. Finally in Sections 



4.5 and |4.6[ we study the mass accretion and luminosity histories of our simulations. 



4.1. General Results 

Figure [l] shows snapshots of the distribution of cloud tracer material for model cc_il_bl_95 over 10 years of its orbit, 
starting in 2010. We do not show earlier epochs, since the cloud is largely spherical prior to about 2010. As we are 
presenting three-dimensional simulations, we generally use volume visualizations, which allows one to see the entire 
simulation domain. There is, of course, an opacity factor, so one can only see through low tracer regions. Only in 
r igures [2] and J3] do we show two-dimensional slices through our data. 

As expectedirom our estimates in Section |3) the early phase of evolution is dominated by tidal stretching. By 
2012.5, the cloud covers an arc of roughly 1.0 x 10^^ cm (5.3i?*), compared to widths and depths of roughly 2.0 x 10^^ 
cm (l.li?*) and 1.4 x 10^^ cm (0.73i?*), respectively. Ram pressure forces become important at later phases when they 
compress the head of the cloud and slow its orbital motion. Hydrodynamic instabilities acting throughout the cloud 
history produce streams of material that drift away from the main cloud body, particularly along its inner edge. As 
this material loses angular momentum relative to its Keplerian value (through interaction with the background gas), 
it accelerates toward the black hole ahead of the main cloud body. By 2013.5, this material begins to penetrate within 
the accretion radius. 

Beginning around 2013.5, the destruction of the cloud greatly accelerates, so that, by 2015, it has separated into 
numerous filaments, with additional cloud material dispersed over a volume many times larger than the initial volume 
of the cloud. By 2020 it is clear that very little, if any, cloud material remains on the original Keplerian trajectory. 
Instead, the bulk of the cloud material has been stretched out into an extended filament, anchored at the accretion 
radius race- 

The overall morphology of model cc_il_bl_95 shows some similarities to the comparable two-dimensional model, 
CCOl, presented in Schartmann e t ar| ( |2QT2 ) (cf. their Figure 4 vs. our FigureJT]). Certainly the position of the bulk of 
the cloud material is similar in each corresponding epoch. However, many differences are also apparent, especially in 
the fine- level detail of the accretion filaments. For these, the two-dimensional simulations appear to show a very thin, 
delicate stream connecting the bulk of the cloud to the accretion radius, whereas the three-dimensional simulations 
show a much more dispersed distribution. There are two factors contributing to this difference. One is simply a 
projection effect. If we were to take a single slice through one of our three dimensional simulation volumes, as we do 
in Figure |2j we would see a much greater level of fine detail, which is lost when looking through the entire volume 
with opacity effects. The other, more crucial, factor is that two dimensional simulations severely restrict the avenues 
by which the cloud can break up. The extra degree of freedom in three dimensional simulations has a significant effect 



on the evolution of the cloud, a point we return to in Section 4.5 



4.2. Equation of State Comparison 

Figure [2] shows the temperature of the cloud for our three different EOS models on date 2013.5. The isothermal 
(r = 1) cloud model, cc_il_bl_95, of course maintains its initial temperature T = 10^ K. The isentropic (F = 5/3) 
model, cc_i53_bl_95, shows some temperature variation, but the bulk of the cloud material remains below 10^ K. The 
polytropic (F = 5/3) model, cc_p53_bl_95, on the other hand, shows considerable heating, with much of the cloud 
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Figure 1. Three-dimensional, volume visualization of cloud tracer material T for model cc_il_bl_95, spanning the period 2010 to 2020. 
The small gray sphere appearing in each panel beginning in 2012.5 indicates the accretion volume, r < race = 750Rs. The color represents 
the amount of tracer (value of T) in a particular region, with the transparency also dependent upon the value of T (as indicated by the 
slanted gray bar above the color bar). The cloud material initially has a uniform value of To = 2.36. 
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Figure 2. Pseudocolor plots of temperatures for models cc_il_bl_^5, cc_i53_bl_95, and cc_p^_bi_95 on date 2013.5. Tnese plots are 
two-dimensional slices in the z = plane. The black circle indicates the accretion radius, race = 750Rs, while the white contour indicates 
a cloud tracer level of T = 0.5. 

approaching 10^ K. The differences between the two F = 5/3 simulations indicate that non-adiabatic processes (i.e. 
shocks) play a significant role. This can be seen more directly by noting in Figure p] that the isentropic constant k, 
retains its initial value in models cc_il_bl_95 and cc_i53_bl_95, but changes quite substantially for the more general 
polytropic model cc_p53_bl_95. 
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Figure 3. Pseudocolor plots of the isentropic constant n = P/p^ , normalized to its initial value, for models cc_il_bl_95, cc_i53_bl_95, and 
cc_p53_bl_95 on date 2013.5. These plots are two-dimensional slices in the z = plane. The black circle and white contours have the same 
meaning as in Fig. [2] 

There are also some interesting differences in the morphology of the clouds in each of these simulations. For example, 
Figure |2] shows a progression towards greater cloud disruption by the Kelvin-Helmoltz instability as internal heating 
increases. This is attributed to the smaller density differential between cloud and background gas as the cloud volume 
increases and its density decreases with internal heating. A smaller density differential decreases the timescale for 
Kelvin-Helmholtz, as shown in equation ( p!2|) . However, the total mass lost from the cloud by this mechanism by 
2013.5 is less than 3% for all three EOS models (based on measuring the amount of cloud material no longer contained 
within the contours shown in Figures [2] and [3| . 

The most significant difference exhibited by the three equation of state models is the disruption of the head of the 
cloud as it passes pericenter and begins to decelerate. Ram pressure becomes increasingly important at this time as 
the head of the cloud slows and the tail piles up behind. Ray leigh- Taylor instabilities are also triggered at the head 
of the cloud at this time as the background gas effectively accelerates into the denser cloud material. The timescales 
for Ray leigh- Taylor are shortest for greater acceleration or pressure differential between cloud and background gas. 
Hence, as shown in Figure [2j the polytropic model is least affected by Ray leigh- Taylor disruption due to its ability to 
achieve pressure equilibrium with the background gas. 

The polytropic equation of state model indicates that internal dynamic heating of the cloud is considerable. Nev- 
ertheless, the extreme efficiency of radiative cooling is more than enough to keep the temperature of the actual cloud 
close to 10^ K. The isothermal model should, therefore, represent the most realistic simulation of the cloud's evolution. 

4.3. Background Comparison 

The convectively stable background profile given by equations - Q does not appear to lead to a dramatically 
different cloud morphology at early (2012.5 and earlier) nor late (2015 and later) times, as seen by comparing model 
cc_il_b2_95 in Figure H] with model cc_il_bl_95 in Figure [ij there is, however, some difference in the distribution of 
cloud material around the accretion sphere in 2013.5. Nevertheless, we find that the average mass accretion rate in 
cc_il_b2_95 is only about 10% lower than cc_il_bl_95, as shown in Table [T] and Figure [tI This is perhaps not surprising 
since we artificially reset both backgrounds to their initial profiles after each hydro cycle. This artificially stabilizes 
the convectively unstable background. 

Still, the small differences between the two background models can give us some estimate of the level of sensitivity of 
our results to realistic variations in the true galactic center environment. Specifically, our work probes the sensitivity 
to density and temperature. However, other properties of the galactic center may also be important. For example, 
radiation and magnetic fields might affect the evolution of G2, and neither is considered in this work. 

4.4. Start Date - 1944.6 vs. 1995.5 

All but one of our simulations began with a start date of 1995.5, shortly before the discovery of G2. For the other, 
we started the cloud from what would be the apocenter of its present orbit. The cloud in this simulation, cc_il_bl_44, 
obviously has much longer to evolve and traverses a much greater distance through the galactic center. In Figure [5) we 
show what the cloud from this model would have looked like on 1995.5. It is clear that tidal forces are already beginning 
to stretch the cloud and ram pressure forces are blunting its leading edge, as expected based on our calculations in 
Section [3j Clearly if G2 started as a spherical cloud at apocenter, then its shape should have changed significantly 
by the time of its discovery. Ram pressure and hydrodynamic instabilities are also stripping some material from the 
edges of the cloud, generating a tail of tracer material. The mass loss to this point, however, is negligible; the cloud 
still retains 99.95% of its material. Nevertheless, the earlier seeding of these instabilities ultimately leads to a greater 
disruption of the cloud, starting at pericenter passage (cf. Figure [5|. 



Figure 4. 



Same as Fig. ^ but for the alternate background model, cc_il_b2_95. 
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Figure 6. Plot of the total amount of cloud material accreted through race as a function of time, represented as a fraction of the initial 
cloud mass, for four models. 



Our results support the conclusions of Burkert et al. (2012) and "S chart mann et al. (2012) that G2 must have either 
formed shortly before its discovery or be part ot a larger, more extended structure. Otherwise, it should have had a 
much more elongated appearance at the time of its 2011 observation. It is also clear that G2 will not survive even a 
single passage through the galactic center; the environment is just too hostile for a cloud of its size. 

4.5. Mass Accretion 

During the simulations, we trace the mass inflow through the accretion radius, race = 7bORs- Only the accreted 
mass that was originally part of the cloud is recorded. The resulting total accreted mass is plotted against time in 
Figure [g] as a fraction of the initial cloud mass for four of the models (we leave off the isentropic model to prevent 
overcrowding). Cloud material first begins to accrete in 2013. By 2020, roughly 4% of the cloud material has accreted 
in most simulations, giving a baseline mass accretion rate of 5 x 10~^ Mq yr~^. 

Time-resolved mass accretion rates are shown in Figure [7j Averages and standard deviations are also reported in 
Table [l] During the simulations, the time-averaged accretion rate exhibits only stochastic changes over the seven 
year period we monitor. The mean accretion rate in all models re mains fairly constant, and for model cc_il_bl_95 is 



comparable to the analogous two-dimensional model, CCOl, from Schartmann et al. (2012). However, whereas the 



two-dimensional simulations show variability up to an order of magnitude on timescales as short as ~ 1 month, we 
observe variations no larger than about a factor of three on similar timescales. This is another example of where the 
detailed evolution differs significantly in two versus three dimensions. 

One goal of this paper is to estimate how this mass accretion will affect the future activity of Sgr A*. To do this, 
we need to compare our measured accretion rates against the curr ent estimated mass accretion rate in the galactic 
center. Studies (e.g. Blandford fc Begelman|p"999 Yuan et "ar]|2003 ) suggest the following radial dependence 



M{r) = MBondi 



^Bondi 



0.27 



(18) 



with MBondi = 10~^ Mq yr~^ (Yuan et al. 2003) at reondi- This implies that the current mass accretion rate at r = race, 
where we make our measurenients, should be 3.8 x 10~^Mq yr~^. In this case, the capture of G2 material onto Sgr 
A* wo uld only boo s t the observed mass accretion rate by about 1-5%, depending on which of our models we consider. 
The Yuan et al. (2003) estimate, though, assumes that the gas accretes as a radiatively inefficient accretion flow 
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Figure 7. Plot of the instantaneous mass accretion rate through race for four models. The data are sampled every 100th cycle of the 
hydrodynamics solver, giving an effective time sampling of ~ 0.03 yr = 10 d. 

(RIAF). It could be that the cloud material, being cooler, denser, and having a coherent angular momentum, may 
accrete much more efficiently onto the black hole. If that is the case, then our observed mass accretion rate, 5 — 19 x 
10~^Mq yr~^, may lead to a significant increase in the mass accretion rate onto Sgr A*, currently estimated to be 
in the range 2 x 10"^ < M < 2 x 10"'^ yr-^ (Aitken et aL 2000; Bower et aL„2Q03, ,Mar rone et al. 2007). Since 
modeling the galactic center mass accretion rate relies on a number of uncertam assumptions about energy transport 
and radiative efficiency, it will be extremely useful to monitor the break up of G2 to better understand how Sgr A* is 
fed. 

4.6. Luminosity 

We can get a crude estimate (upper limit) of the luminosity associated with the disruption of G2 by tracking the 
amount of cloud energy (gravitational potential plus kinetic) that is dissipated during our simulations. Assuming 
that all of this energy is radiated away instantaneously, we can track its luminosity as a function of time as shown in 
Figure [s] for mo del cc_il_bl_95. The re sulting value, ~ 10^^ erg s~^, is roughly consistent with the predicted pericenter 



luminosity from Gillessen et al. (2012) and should be easily observable. 

5. CONCLUSIONS 

The cloud G2 is currently approaching Sgr A* and in 2013.5, will pass it at a pericenter distance of 3100 times the 
Schwarzschild radius, corresponding to 4.0 x 10^^ cm. This provides a unique opportunity to investigate directly the 
disruption of a cold gas clump by its gravitational interaction with the black hole and hydrodynamic interaction with 
the galactic center e nvironment. Clear evidence for tidal velocity shearing and stretching has already been detected 
dBurkert et"ar]|2012D . 



In this paper we presented five different three-dimensional, numerical, moving mesh simulations illustrating possible 
future evolutionary tracks for G2. Our results generally support previous conclusions (Burkert et al.|20l"2 Schartm ann 



[et al. 201 2) that G2 must have been formed recently or must be only a small piece of a larger cloud conaplex. In any 
case, a cloud like the ones we simulated would not be expected to survive long (apparently not even a single pericenter 
passage) within the galactic center. 

In all of our simulations we find a roughly similar morphology for the cloud: Tidal stretching and Kelvin-Helmholtz 
effects become evident first; then, around the time of pericenter passage, ram pressure and Rayleigh-Taylor instabilities 
additionally help to disrupt the cloud, which begins to fragment; as the cloud disperses over a larger volume, it begins 
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Figure 8. Light curve for model cc_il_bl_95, estimated from the dissipation of the gravitational potential plus kinetic energy of the cloud. 



to lose significant amounts of angular momentum through its interactions with the background gas; the cloud stretches 
into a stream of gas feeding the black hole (for us the accretion volume r < race) at 5 — 19 x lQ~^Mr7^ Thi s 



is consistent with the values seen in the comparable two-dimensional simulation CCOl of Schartmann et al. (2012) 
although the variability of that simulation was up to a factor of three larger. Looking in more detail, we note 
some important differences among our models. The non-isothermal models, particularly the polytropic one, showed 
significantly higher cloud temperatures than the more realistic isothermal ones (which accounted for radiative cooling 
of the cloud). Furthermore, the isentropic model yielded a significantly higher mass accretion rate than the isothermal 
(and even polytropic) models. The model with the denser (convectively stable) background gave a somewhat smaller 
mass accretion rate than the default background. Finally, not surprisingly, model cc_il_bl_44 illustrated that starting 
the simulation earlier leads to greater disruption of the cloud by the time of pericenter passage. 

It is uncertain how much the disruption of G2 will affect the mass accretion feeding Sgr A*. The accretion rate 
observed in our simulations is comparable to the currently estimated rate in the immediate vicinity of Sgr A*. But 
at the radius where we are actually measuring the accretion, it only represents an increase of 1-5% over the currently 
estimated fee ding rate. Such a small cha nge in M may not produce any measurable changes in the Sgr A* emission 
(although see [Moscibrodzka et al. (2012) for examples of what larger accretion rate changes might do to the X-ray, 
infrared, and millimeter emission from Sgr A*). Ultimately it will depend on how efficiently the G2 material can 
get from the radius where our simulations cut off down to Sgr A*. Regardless, the break up of G2 will provide 
an unprecedented opportunity to study accretion physics in the galactic center. Furthermore, our crude luminosity 
estimate of ~ 10^^ erg s~^ suggests this event should be easily observable over the next several years. 

Concerning possible future improvemen ts to our work, we have not ex plored the following alternate scenarios: the 
^^spherical shell" of iBurkert et al. (2012) and Schartmann et al. (2012), the ^^h idden-star /protq - plane tary disk" of 
Murray-Clay fc hoe^^T^ nor the ''photoevaporating, disturbed stellar disk" of iMiralda-Escudel (|2012|) . In the first, 
G2 is just the tip ot a much larger, more massive shell of gas. This scenario may be required to explain the lower 
surface brightness, cone-like structure that accompanies G2 (Gillessen et al. 2012). The second scenario proposes that 
G2 is a dense, proto-planetary disk bound to a low-mass star, which was scattered from the disk of young stars orbiting 
Sgr A*. The third suggests instead that G2 is a cloud of photoevaporated gas originating from a disk around a star 
that was previously disrupted by a stellar mass black hole and now produces a cloud at every periastron passage. 
We have also not considered thermal conduction. This is likely the next most important physical influence beyond 
those considered in this work. Thermal conduction will be especially important at late times as the cloud breaks up. 
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Magnetic fields may also play an important role. For exan iple, magnetic fields can affect the growth and development 
of hydrodynamic instabilities along the cloud surface (e.g. Fragile et al.||2QQ5" ). 
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APPENDIX 
MOVING MESH 

As this is the first introduction of moving meshes implemented in the Cosmos++ code, we present a brief description 
of the covariant framework from which the equations are derived. A more formal presentation of numerical methods 
and code tests will be considered in another paper. Throughout this Appendix we use standard index notation 
in which repeated indices represent summations over spatial components, and the raising and lowering of indices (to 
contravariant and covariant components) is done with the spatial metric, e.g. Vcx = Qa^v^ - We also adopt the convention 
in non-relativistic work of using Latin indices {i^j^k) to represent quantities measured in Cartesian coordinates, and 
Greek indices (a,/3,7) for generalized curvilinear coordinates. The indices run over the three spatial dimensions (they 
do not include the time component). 

Before presenting the equations, we draw an analogy with general relativistic magnetohydrodynamics (MHD) to 
define state variables and write the stress-energy tensor for a perfect fluid (and ideal MHD) as a linear combination 
of the hydrodynamic T^^ and magnetic contributions T^^: 

j.a(3 ^ j.a(3 ^ j.a(3 ^ ^^a^(3 ^ (p ^ Pb)^/"^ + Q"^ - b^'b^ , (Al) 

where p is the fluid mass density, is the contravariant fluid velocity, P is the fluid pressure (for an ideal gas 
P = (r — l)e where e is the fluid internal energy density and F is the adiabatic index), is the magnetic field, 
Pb = ^a/3^"^^/2 is the magnetic pressure, is the tensor artificial viscosity used for capturing shocks, and g^^ 
is the spatial 3-metric tensor for general curvilinear coordinates. Although we do not model magnetic fields in this 
paper, we include the full MHD equations for completeness and future reference. 

We begin by writing in flux- conserving form the evolution equations for mass, internal energy, momentum, and 
magnetic induction in Cartesian coordinates (x*) as: 

dt + dx^ dx^ dx- ' ^ 

where 5* is the Kronecker delta tensor, and = gij is effectively the flat metric in Cartesian coordinates. Additionally 
is a scalar potential introduced as a divergence cleanser to maintain a divergence- free magnetic field (pib^ = 0), and 
(\) is the gravitational potential satisfying Poisson's equation 

(|^ = 4^G.- (A6) 

Next we consider an arbitrary time-dependent transformation into generalized coordinates with x(^,t), triad 
= spatial 3-metric go,^ = e^e^^ij, metric determinant = Det(| |^q,/3| |), and grid velocity i:* = dx'^ /dt\^ = 

Vg. Applying standard transformation rules for gradients and vector quantities in general, the MHD equations take 
the following form in moving generalized coordinates: 

^ + l^TspenV - V;)] = , (A7) 

^ + [vs-f (..• - K)] = - + 03) ^ . (A8) 



Ot • g^a - V^)] =-— {yTge^ [{P + P,)51 + Q1- 6-6,] } 



V^PS^^, (A9) 
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dt 



[V9t>>ef{v'-V;)]=^b^ 



V99' 



(AlO) 



(All) 



In this form the evolved vector fields (velocity, momentum and magnetic field) are defined in Cartesian coordinates 

Taking these transformations a step further, we assume the triad matrix is a product of two transformation bases 

(A12) 



where is a static (time independent) transformation between Cartesian coordinates and any other fixed frame 

coordinate set (^"), for example spherical or cylindrical coordinates. The second triad is a dynamical (time 
dependent) transformation between the static fixed frame coordinate system (spherical, cylindrical) and generalized 
curvilinear coordinates. This decomposition allows for grid motion in more general topological "background" grids 
besides Cartesian. It affects only the vector conservation equations, allowing vector quantities to be e xpre ssed i n the 
most convenient (symmetry respecting) triad. Hence only the momentum and magnetic field equations (A9) and (AlO) 
are modified under this transformation as 



'Vds^e^v^ - Vp] =-g- { v^e| [(P + P,)5^ + - b^b,] } 



dt 



(A13) 
(AM) 



The mixed index (hatted and unhatted) term T^^^ in equation (A13) is actually Sije'^dae^-^^ and T^^ in (|A14|) is 



e^^c^e^:, which follow from the Christoffel/triad gradient identities. Both of these terms are time-independent and can 
be evaluated as static triad Christoffel symbols at the initial time when dynamical and static triad coordinate systems 
are identical. 
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